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The coupling of slosh dynamics within a partially filled rotating dewar of superfluid 
helium II with spacecraft orbital dynamics is investigated in response to the environ- 
mental disturbances of (a) lateral impulses, (b) gravity gradients and (c) g-jitter forces. 
The purpose of this study is to investigate how the coupling of helium II fluid slosh 
dynamics driven by three cases of environmental force with spacecraft dynamics can 
affect the bubble deformations and their associated fluid and spacecraft mass centre 
fluctuations. The numerical computation of slosh dynamics is based on a rotational 
frame, while the spacecraft dynamics is associated with a non-rotational frame. Results 
show that the major contribution of orbital dynamics is driven by coupling with slosh 
dynamics. Neglecting the effect of slosh dynamics acting on the spacecraft may lead 
to the wrong results for the development of orbital and attitude control techniques. 
© 1996 Elsevier Science Limited 
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In order to carry out scientific experiments, some exper- 
imental spacecraft use cryogenic cooling for telescope 
instrumentation, and occasionally to maintain temperatures 
near absolute zero for mechanical stability. Most fre- 
quently, the cryogen used is helium II. In this study, the 
coupling of spacecraft orbital dynamics with sloshing of 
helium II driven by lateral impulses, gravity gradients and 
g-jitter forces is investigated. Both gravity gradients and g- 
jitter forces will be defined later. In the specific case of the 
Gravity Probe-B (GP-B), the spacecraft uses liquid helium 
as a cryogen and propellant to maintain the cooling of 
instruments, attitude control and drag-free operation of the 
spacecraft. Potential fluid management problems may arise 
due to an asymmetric distribution of liquid helium and 
helium vapour due to perturbations in the liquid-vapour 
interface. A basic understanding of the coupling of slosh 
dynamics in six degrees of freedom with orbital dynamics 
plays a significant role in the development of spacecraft 
guidance and attitude control systems. 

Liquid helium at a temperature of 1.8 K is used on the 
GP-B. Because it is a superfluid, there are at most only 
very small temperature gradients in the liquid helium. In 
the absence of Marangoni convection, due to the negligibly 
small temperature dependence of surface tension and negli- 
gible thermal gradients along the liquid-vapour interface, 
the equilibrium shape of the interface is governed by a bal- 


ance of capillary, centrifugal, gravitational and dynamic 
forces. Determination of liquid-vapour interface profiles 
based on computational experiments can uncover details of 
the flow that cannot be easily visualized or measured exper- 
imentally in a microgravity environment. 

An instability of the liquid-vapour interface can be 
induced by the presence of longitudinal and lateral acceler- 
ations. Thus, slosh waves are excited, producing high- and 
low-frequency oscillations in the liquid helium. The sources 
of the residual accelerations include effects of the Earth’s 
gravity gradient and g-jitter 13 . A recent study 4 suggests 
that the high-frequency accelerations may be unimportant 
in comparison with the residual motions caused by low- 
frequency accelerations. 

The time-dependent behaviour of liquid helium in a par- 
tially filled rotating dewar in reduced gravity environments 
was simulated by numerically solving the Navier-Stokes 
equations subject to initial and boundary conditions 5,6 . At 
the interface between the liquid and vapour, both the kine- 
matic surface boundary condition and the interface stress 
conditions for components tangential and normal to the 
interface were applied 2 ' 3 . The initial conditions were 
adopted from the steady-state formulations developed by 
Hung et al . 1 . Some of the steady-state formulations of inter- 
face shapes were compared with the available experiments 
carried out by Leslie 8 in a free-falling aircraft (KC-135). 
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The experiments of Mason et a/. 9 showed that classical 
fluid mechanics theory is applicable to cryogenic helium in 
large containers with sufficiently large velocities 10-13 . 

At temperatures close to absolute zero, quantum effects 
begin to be of importance in the properties of fluids. At a 
temperature of 2.17 K, liquid helium undergoes a second- 
order phase transition; at temperatures below this point, 
liquid helium (helium II) has a number of remarkable 
properties, the most important of which is superfluidity. 
This is the property of being able to flow without viscosity 
in narrow capillaries or gaps. At temperatures other than 
zero, helium II behaves as if it were a mixture of two differ- 
ent liquids. One of these is a superfluid and moves with 
zero viscosity along a solid surface. The other is a normal 
viscous fluid. The two motions occur without any transfer 
of momentum from one to the other for velocities below a 
critical velocity 10 12 . Above the critical velocity, the two 
fluids are coupled through their mutual friction 10 L \ 

The key parameters affecting critical velocity are tem- 
perature and container size. To determine the dynamic 
behaviour of helium II in a large rotating cylinder, the mut- 
ual friction of the two fluid components is accounted for 
in the model computation. The density concentration of the 
superfluid is a function of temperature, which is also true 
for the surface tension and viscous coefficient for helium 
II 10 17 . In this study, the theory of viscous Newtonian fluids 
is employed with transport coefficients being a function 
of temperature. 

In order to carry out the study of transient phenomena 
of coupling between slosh reaction torques driven by three 
types of acceleration and spacecraft orbital dynamics, (a) 
slosh dynamics based on a fluid dynamics formulation and 
(b) orbital dynamics based on translational and rotational 
formulations of spacecraft dynamics, have been numeri- 
cally solved simultaneously. The three types are (a) lateral 
impulse, (b) gravity gradient-dominated and/or (c) g-jitter- 
dominated accelerations. In other words, coupling between 
slosh dynamics and orbital dynamics will deeply affect the 
behaviour of the bubble driven by slosh dynamics, which 
is quite different from the results obtained earlier 15,18 with- 
out the implementation of the orbital dynamics modifi- 
cation. Meanwhile, it is true that orbital dynamics driven 
by slosh dynamics can cause the spacecraft to deviate from 
normal operation 19 . 


Non-inertial frame mathematical formulation 
of slosh dynamics 

An experiment by Andronikashvili 10-13 for rotating helium 
II showed that it is necessary to exceed a critical velocity 
for the interaction between the normal and superfluid 
components to establish rotation of the entire bucket 10-13 . 
For a rotating dewar with outer diameter 1 .56 m and inner 
diameter 0.276 m, the critical velocities are 6.4 x 10 -7 and 
3.6 x I0 -6 m s -1 , respectively 10-15 . With a rotating speed of 
0.1 rev min -1 , the linear velocities along the outer and inner 
walls of a rotating dewar are 8.17xl0 -3 and 
1.45 x 10 -3 ms -1 , respectively, which are at least several 
hundred times greater than the corresponding critical velo- 
cities. Based on this illustration, the problem under con- 
sideration has the necessary features for a viscous Newton- 
ian fluid formulation to be adopted in this study. 

Consider a closed circular dewar partially filled with 
helium II. The whole fluid system is spinning in the axial 


direction z of cylindrical coordinates (r, 6 , z), with corre- 
sponding velocity components (m, v, vv). The governing 
equations for non-inertial frame bound coordinates of a 
spacecraft spinning along its z-axis have been given in 
recent studies 20-22 . Thus, dynamic forces relating to gravity 
gradient, g-jitter and angular accelerations, and centrifugal, 
Coriolis, surface tension, viscous forces, etc., are given 
explicitly in the mathematical formulations 18 22 . In the 
computation of sloshing reaction forces, moments, viscous 
stress and angular momentum acting on the container wall 
of the spacecraft, one must consider those forces and 
moments in the inertial frame rather than in the non-inertial 
frame 20 24 . 

For the purpose of solving slosh dynamic problems of 
liquid systems in orbital spacecraft under a microgravity 
environment, one must solve the governing equations 21 - 23 
accompanied by a set of initial and boundary conditions. 
A precise and detailed illustration of these initial and 
boundary conditions was given by Hung and Pan 19-22 24 . 
The computational algorithm applicable to cryogenic fluid 
management under microgravity was also given earlier 20 
24 . Summarized computational algorithms are illustrated in 
Figure 1 . In this study, in order to show a realistic example, 
a dewar with an outer radius of 0.78 m, and an inner radius 
of 0.138 m, top and bottom radii of 1.10 m and a height of 
1 .62 m has been used in the numerical simulation. The 
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dewar tank is 80% filled with liquid helium and the ullage 
is filled with helium vapour (total fluid mass 287.6 kg). 
The temperature of liquid helium is 1.8 K. In this study the 
following data were used: liquid helium density 

= 145.7 kg m~\ helium vapour density - 0.45 kg m \ fluid 
pressure = 1 .6625 Pa, surface tension coefficient at the 
interface between liquid helium and helium vapour 
= 0.0353 N rrr\ liquid helium viscosity coefficient 
= 9.61 x 10~ 9 m 2 s _1 and contact angle = 0°. The initial pro- 
files of the liquid-vapour interface for the rotating dewar 
are determined from computations based on algorithms 
developed for the steady-state formulation of microgravity 
fluid management 7,23 . 

A staggered grid for the velocity components is used in 
this computer program. The MAC (marker-and-cell) 
method 26 of studying fluid flows along a free surface is 
adopted and the VOF (volume of fluid) method is used to 
solve finite difference equations numerically. The approxi- 
mate flow velocity is calculated from the explicit approxi- 
mation of momentum equations based on the results from 
the previous time step. Computation of pressure and velo- 
city at the new time step are thus obtained from iteratively 
solving the pressure equation through the conjugate 
residual technique 27 30 . The configuration of the liquid- 
vapour interface adjusted by the surface tension effect at 
the new time step is then obtained. The time step during 
this computation is automatically adjusted through the 
fulfilment of the stability criteria of computed grid size. 
The convergence criterion of the iteration of the pressure 
equation is based on the computed velocity at each cell 
which satisfies the continuity equation with errors of no 
more than 10 -3 in the velocity difference 31 . As for the vol- 
ume conservation of liquid, a deviation of less than 1% 
error in volume is guaranteed before moving to the next 
time step. 

In this study, characteristics of the slosh reaction forces 
and torque fluctuations exerted on the dewar in response 
to various accelerations acting on the spacecraft are also 
investigated. The mathematical formulation of the fluctu- 
ations of slosh reaction forces and torques exerted on the 
dewar are illustrated elsewhere 21,32,33 . With reference to 
slosh dynamics driven by impulse 13,19 and mathematical 
formulations 21,32,33 , one can calculate the slosh reaction 
force and the associated torque acting on the dewar, calcu- 
lations that are required for the computation of the coupling 
of the slosh dynamics of He II with the orbital dynamics 
illustrated in the next section. 


Inertial frame mathematical formulation of 
spacecraft dynamics 

In spacecraft dynamics, a rigid body with six degrees of 
freedom, three being translational and three rotational, is 
considered. In this study, our primary interest is to investi- 
gate the coupling between slosh reaction torques driven by 
(a) lateral impulse, (b) gravity gradient-dominated and (c) 
g-jitter-dominated accelerations and spacecraft dynamics. 
In other words, our main purpose of study is to investigate 
how the coupling of slosh reaction forces/torques with 
spacecraft dynamics can affect cryogenic bubble defor- 
mations and their associated fluid and spacecraft mass cen- 
tre fluctuations. 

The three governing translational equations are given by 


£ (!»*,) = F n + F u ( 1 ) 

dr 

where w, X, F D and F { denote the mass of the spacecraft, 
the inertial frame (non-rotational) coordinate, the residual 
environmental force acting on the spacecraft and the slosh 
reaction force (from the fluid system) acting on the space- 
craft, respectively 19,20,30,31 . Subscript / denotes components 
along the /-direction (= X, Y or Z in the non-rotational 
frame), while single and double overdots imply first- and 
second-order time derivatives of a parameter, respectively. 
Eulerian angles are defined to accommodate three rotational 
equations 34 . As usual, i//, 0 and </> are defined as the heading, 
attitude and bank angles 34 . Three rotational equations in 
terms of Eulerian angles are given by 

/ t (0 + iff sin 0 cos0) + l. (i ip cos 0(</) - *//sin0) = M 0 (2) 
/ a [(-</H- ips'md+ijjO cos0)sin0 - 6 K<f>~ </' sin/)) cos0] + 
/t(i//cos 2 0 - 2ifjb sin0 cos0) = (3) 

/ a (</>- i//sin0- i//0 cost)) = M# (4) 

where 


M 0 = (Mi t + M Uv ) cos 4> - (M l v + M Dv ) sine/) (5) 

M (/ , = (M l x + M d J sin</> cos 0 + 

(M Lv + M Dv ) cose/) cos 0- (M u . + M Dr ) sin 0 (6) 

M+ = (M u + M Dz ) ( 7 ) 

Here M u and M t>( are the slosh reaction torque acting on 
the dewar and the residual torque acting on the spacecraft, 
respectively, along the /-direction (/= x, y or z in the non- 
inertial frame). / a (=/,-) and I { ( = / tA =/ vv ) denote the 
moments of inertia along the axial and transverse direc- 
tions, respectively. 

Activation of the coupling of slosh dynamics in response 
to lateral impulse and environmental disturbances, such as 
gravity and g-jitter accelerations, with spacecraft dynamics 
certainly will induce angular velocities (co^a^o**) along 
with (x,y,z) coordinates in the rotational frame. We can 
assume that (to> 1 ,w 2 ,<u 3 ) = (0j,0 2 ,03 + w). Here (0j,0 2 ,0 3 ) are 
the time derivatives of the angular displacements along the 
(x,\\z) coordinates and co denotes the rotational speed along 
the i-axis. With the definitions / VA . = / ,, f YV = l 2 and l :: - / 3 , 
the spinning angular velocity in the rotational frame can be 
expressed as follows: 

e , + k 2 o 2 ( b, + (*)= ] - (M u + m D v ) 

^ i 

02 + x 2 0,(0 3 + 0 >)=J (M Lv + M Dv ) (8) 

l 2 

0 3 + KA 02 — ~r ( 37 l - + M 0z ) 

h 

where X, = (/ 3 - / 2 )//„ K 2 = (/, - / v )// 2 and X 3 = (I 2 - /, )// 3 . 

Both the translational and rotational equations are initial- 
value problems. Initial conditions for the translational and 
rotational equations are 


(X, Y, Z) = (0, 0, 0) cm 
(X, KZ) = ( 0, 0, 0) cm s _l 


at t = 0 


(9) 
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( 0 , i/>, </>) = (0, 0, 0) rad 1 

( 0, if/, (f > ) = (0, 0, 0. 1 ) rev min 1 J at r ^ 1 * ^ 

In order to solve the coupling of slosh dynamics and space- 
craft orbital dynamics, as shown in Equations (1) to (4), 
an iteration method at every time step was adopted to deter- 
mine (F d „ M Di ) for spacecraft that are mutually coupled 
with {F u% M u ) from the slosh dynamics. The summarized 
computational algorithm is illustrated in Figure 2, where 
the superscript k is the number of the iteration. In this com- 
putation, the dry mass of the spacecraft (excluding the fluid 
mass) is assumed to be 350 kg, the fluid mass with an 80% 
liquid-filled level is 287.6 kg and the moment of inertia for 
the spacecraft is /„ = / vv = 160.3 kg m 2 , and /„ = 

154.5 kg nr. 


Coupling of He II slosh dynamics and orbital 
dynamics 

The characteristics of the three types of force on fluid sys- 
tems are quite different. The lateral impulse is in the form 
of a spotty 6-function with a short time period of 10 -2 s, 
while both the gravity gradient and the g-jitter forces are 
in the form of a continuous spectrum covering the whole 
spacecraft orbit period. The magnitudes of these acceler- 
ations vary from 10" 2 to 10“ 9 g o (where g 0 =9.81 ms -2 ) 
and act on fluid elements at different locations either with 
different magnitudes (gravity gradient forces) or the same 
magnitudes (impulse and g-jitter forces), depending on 
their characteristics. These are illustrated in Table J with a 
full description shown in this section. 



Figure 2 Computational algorithm for coupling of slosh and orbital dynamics 
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Table 1 Characteristics of various orbital accelerations driven slosh dynamics 



Lateral impulse 

Gravity-gradient- 

dominated 

G-jitter-dominated 

Types of forces 

Spotty 6-function 

Continuous spectrum 

Continuous spectrum 

Magnitudes of forces 

10 3 g Q 

10 7 g„ 

10- 6 9o 

Time period of force acting 

10 2 s 

Full orbit period 

Full orbit period 

Magnitude ranges of forces 

(10 2 to 10~ 9 ) g Q 

10" 7 g 0 

(10 5 to 10 9 ) g Q 

Magnitudes of force acting on fluid element at 
different locations 

Same magnitude 

Different magnitudes 

Same magnitude 

Types of fluid motions 

Leftward and 
rightward oscillations 

Tidal mode motion 

Leftward- rightward 
and up-down 
oscillations 

Maximum fluid mass centre fluctuations (cm) 

2 

1.4 

5 


Coupling between lateral impulse driven slosh 
dynamics and spacecraft dynamics 

It is assumed that a lateral impulse with the following forms 
of force (F,) and torque (A/,) acts on the spacecraft system 
in non-inertial (rotating) frame Cartesian coordinates 1518 : 


F, = (F x ,F y ,F z ) = ( 20,0,0)N 1 

M, = (A/,, M y , M z ) = (0, 10, 0) N mj 


for 0 < t < 1CT 2 s 

(ID 


and 


F i — ( 0, 0, 0) N 
Mi = (0, 0, 0) N m 


for t > 1 0 2 s 


( 12 ) 


For a spacecraft mass of 637.6 kg, the lateral impulse is 
equivalent to 3.2 x 10 ~ y g 0 in this case. In general, the mag- 
nitude range of the lateral impulses is I0~ 2 to 1 0“ 9 g (V 

The entire computation of the slosh dynamics is carried 
out for coupling with the spacecraft dynamics. Slosh 
dynamics concluded from the present study are quite differ- 
ent from the results obtained earlier 15 - 18 without the con- 
sideration of spacecraft dynamics. Figure 3 shows the time- 
sequence evolution of the three-dimensional dynamic 
behaviour of the liquid-vapour interface (bubble) oscil- 
lations (in the rotational frame) activated by a lateral 
impulse acting on the spacecraft in the positive ^-direction 
and reacting to the fluid in the negative jc-direction (liquid 
in negative x -, and bubble in positive x-directions). The 
time-sequence evolution illustrates the following behaviour 
of the bubble dynamics in response to a lateral impulse: 
(a) the bubble was first shifted to the positive jc-direction 
(liquid to the negative jc-direction); (b) the bubble was also 
shifted to the positive y-direction due to the Coriolis force 
with rotation along the z - axis in the rotational frame; and 
(c) bubble deformations with back-and-forth oscillations 
are viewed during the 1500 s simulation time span. 

Figure 4A shows the time evolution of the growth and 
decay of fluid mass centre fluctuations (in the rotational 
frame) in response to a lateral impulse. It is also due to the 
torque of the impulse acting on the y-axis which drives the 
propagating slosh dynamic disturbances to the z-axis. The 
values of the maximum and minimum locations of the fluid 
mass centre fluctuations are Max (jc lo y LC , Zlc) = (0.46, 
10 -4 , 0.25) cm and Min (* LC , y LC , z LC ) = (-1.3, -1.9, 
-0.28) cm, respectively, while the fluctuating values of the 
maximum fluid mass centre disturbances are Max (Ax LC , 
Ay LC , Azlc) = (1.73, 1.93, 0.53) cm. This shows that Ay LC 
> Ajc lc > Az lc . The response to a lateral impulse can be 


described as follows: (a) the fluid mass was shifted to the 
negative jc-direction initially due to the dewar reaction force 
acting on the fluid in response to the impulse on the space- 
craft in the positive .^-direction; (b) the shift in the fluid 
mass centre fluctuations to the negative y-direction in the 
rotational frame was also caused by the dewar reactional 
Coriolis force acting on the fluid in response to the impulse; 
(c) this shows that the major fluid mass centre fluctuations 
are in the negative jc, negative y quadrant with small ripple 
oscillations in the z-axis; (d) there is a time lag of 150 s 
between the impulse acting on the spacecraft and the 
maximum fluid mass centre displacement; in other words, 
the fluid mass centre reacted rather slowly; (e) Figure 4B 
shows an exponential decay in the time history of the fluid 
mass centre fluctuations with a very long time constant 
because of the extremely low surface tension and viscosity 
of helium II fluids. 

The spacecraft mass centre fluctuations (in a non- 
relational frame) were also studied, as illustrated in 
Figure 4C . The maximum and minimum locations of the 
spacecraft mass centre fluctuations are Max (X sc , K sc , Z sc ) 
= (0.85, 0.63, -0.24) cm and Min (X sc , T sc , Z sc ) = (-0.6, 
-0.76, -0.50) cm, respectively, while the fluctuating values 
of the maximum spacecraft mass centre disturbances are 
Max (AX SC , AK SC , AZ SC ) = (1.46, 1.39, 0.26) cm. This shows 
that AX SC > AK SC > A Z sc in a non-rotational frame. The 
spacecraft mass centre fluctuations are very similar to those 
of the fluid obtained in Figures 4A and 4B except that the 
former is in a non-rotational frame, and the latter in a 
rotational frame in which the Coriolis force transfers a dis- 
turbance from the X- to the T-axis. The magnitude is 
smaller by the ratio of the liquid weight to the total space- 
craft weight. Figure 4D shows the time history of the 
spacecraft mass centre displacement. Clearly the effects of 
helium slosh must be considered in the spacecraft attitude 
control systems. 


Coupling between slosh dynamics driven by 
combined gravity gradient and g-jitter 
accelerations dominated by gravity gradient 
forces and spacecraft dynamics 

To give an example of the GP-B spacecraft, which is an 
Earth satellite orbiting at 650 km directly over the poles, 
the orbit period, r () , can be computed from the following 
expression: 


t 0 = 27 T 


/? C V2 


(13) 


where R E denotes the radius of the Earth (=6373 km), R c 
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Figure 3 Three-dimensional bubble deformations in response to coupling between lateral impulse and orbital dynamics 


the radius of the circular orbit (= R r . + h = 7023 km), h the 
orbit altitude (= 650 km) and g 0 the Earth’s gravity acceler- 
ation. For the case of GP-B, the orbit period t 0 = 97.6 min 
(5856 s), and the orbit rate n = 27r/r 0 = 1.07 x 10 -3 rad s 1 . 

The gravity gradient acceleration acting on the fluid mass 
of the spacecraft can be given by 18-20 

& gg = -n 2 [3(r c ’d)r c - dj (14) 

where a gg denotes the gravity gradient acceleration vector, 
d the vector (not a unit vector) from the fluid element to 
the spacecraft mass centre, f c a unit vector in the direction 
from the spacecraft mass centre to the centre of the Earth 
and n the orbit rate. 

In the case of the GP-B spacecraft, it is assumed that the 
gravity gradient exerted on the mass centre (balanced by 


the centrifugal forces of the orbiting spacecraft) is zero 
when the spacecraft is moving around the Earth on its 
specified orbit. In other words, the gravity acceleration 
exerted on the spacecraft consists only of the gravity gradi- 
ent acceleration which is defined in Equation (14). 
Figure 5A illustrates the geometrical relationship between 
the parameters shown in Equation (14). 

At time t = 0, the rotating axis of the spacecraft is aligned 
with the radial direction from the centre of the Earth to 
the spacecraft mass centre. The azimuth angle of the Earth 
towards the location of the spacecraft mass centre, can 
be computed from the orbit period obtained from Equation 
(13) under normal operation of the spacecraft: 

*Ae • — 2 t 

T o 


( 15 ) 
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(A) Fluid Mass Center (Rotational Frame) (B) Fluid Mass Center (Absolute Value) 
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(C) Spacecraft Mass Center (D) Spacecraft Mass Center (Absolute Value) 

(Non - Rotational Frame) 
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Figure 4 Time evolution of fluid and spacecraft mass centre fluctuations driven by coupling between lateral impulse and orbital 
dynamics. (A) Fluid mass centre in (x, y, z) components {rotational frame), (B) fluid mass centre (absolute value), (C) spacecraft 
mass centre in (X, Y, Z) components (non-rotational frame) and (D) spacecraft mass centre (absolute value) 


(A) Geometry of Gravity (B) Time Variation of Gravity (C) Time Variation of Gravity 

Gradient Acceleration Gradient Acceleration Jitter Acceleration 



Figure 5 Time evolution of gravity gradient and g-jitter accelerations acting on the spacecraft during a full orbit period. (A) Geometry 
of gravity gradient acceleration, (B) time evolution of gravity gradient acceleration acting on fluid element located at ( r , 6, z ) = (40 cm, 
7 t /4, 10 cm) and (C) time evolution of g-jitter acceleration 
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where t is the time measured from the instant when the 
direction of the spacecraft’s spinning axis is aligned with 
the radial direction from the spacecraft mass centre to the 
centre of the Earth. Some modification is required if the 
slosh and orbital dynamics are coupled and the results for 
the spacecraft orbit deviating from normal operation are 
considered. 

Fluctuations in the residual gravity due to g-jitter acceler- 
ation are modelled by the following equation 18-20 : 


£ = £b 


1 + ^ sin(2vft) 


(16) 


where g B denotes the background gravity environment and 
/(Hz) stands for the frequency of the g-jitter. 

In this study the following values have been considered 4 
for the investigation of the oscillations of the liquid-vapour 
interface: g-jitter acceleration with a background gravity of 
10 -8 g 0 due to spacecraft atmospheric drag, a rotating speed 
of 0.1 rev min -1 for normal GP-B spacecraft operation and 
0.1, 1.0 and 10 Hz for the g-jitter frequencies. Components 
of the g-jitter acceleration in the non-inertial coordinate 
system are given by 18-20 

_ 

= [g sin ijj E cos (0 + a>/), - g sint \jf v sin ( 0 + cot) 

-g cost/^] (17) 

The characteristics of the gravity gradient and jitter acceler- 
ations, shown in Equations ( 14) and (17), respectively, are 
quite different. The gravity gradient acceleration indicates 
the following two characteristics: (a) the acceleration acting 
on any fluid mass inside the container increases by two 
units of acceleration per unit of distance measured from the 
mass centre of the container (point O g in Figure 5A) to the 
location of the fluid mass parallel to the radial axis from 
the mass centre of the container to the centre of the Earth 
(parallel to the unit vector r c shown in Figure 5A); and (b) 
the acceleration acting on the fluid mass decreases by one 
unit of acceleration per unit of the shortest distance meas- 
ured from the location of the fluid mass to the radial axis 
along the vector from the mass centre of the container to 
the centre of the Earth 17-19 . As an example, Figure 5B 
shows the time variation of the gravity gradient acceler- 
ations for a full orbit period of 5856 s with container rotat- 
ing speed of 0.1 rev min -1 for components along the (x, 
v, z) directions acting on the fluid mass located at (r, 0, 
z) = (40cm, 7t/4, 10 cm). As the magnitude and direction 
of the gravity gradient acceleration acting on each fluid 
mass is strongly dependent on how far the location of the 
fluid mass deviates from the mass centre of the container 
measured along the axis parallel to the vector f C5 which 
varies with respect to time, it shows that the gravity gradi- 
ent acceleration acting on the fluid mass is different for 
fluid masses at different locations in the container. The 
characteristics of symmetry are shown between the first and 
second half-orbit periods of the time evolution of gravity 
gradient acceleration if the orbit period is on an integral 
number of spacecraft spinning periods and the spacecraft is 
under normal operation 18 20 . Coupling of slosh and orbital 
dynamics resulted in Figure 5B and shows that the time 
evolution of the gravity gradient acceleration is asymmet- 
ric, whereas it is symmetric for the non-coupling case 18 20 


Figure 5 B also shows that the magnitude of the gravity 
gradient accelerations is of the order of 10 -7 g o . Unlike 
gravity gradient accelerations of the fluid mass which have 
different magnitudes at different locations in the container, 
the g-jitter acceleration provides the same acceleration to 
the fluid mass at different locations in the container. 
Figure 5C shows the time variation of g-jitter accelerations 
for a full orbit period of 5856 s with a container rotating 
speed of 0.1 rev min -1 and a jitter frequency of 0.1 Hz for 
components along the (jc, y, z)-directions acting on the fluid 
mass everywhere in the container. Again, symmetry is dem- 
onstrated between the first and second half-orbit periods of 
the time evolution of the g-jitter acceleration if the orbit 
interval is an integral multiple of the spacecraft spinning 
period and the spacecraft is under normal operation 18 20 . 
The coupling of slosh and orbital dynamics resulted in 
Figure 5C and shows that the time evolution of the g-jitter 
acceleration is asymmetric, whereas it is symmetric for the 
non-coupling case 18-20 . With the background gravity of 10 -8 
g 0 for g-jitter acceleration considered in this study, the com- 
bined gravity gradient and g-jitter accelerations acting on 
the spacecraft correspond to the case of orbital acceler- 
ations dominated by gravity gradient forces (of the order 
of 10 -7 g G ) acting on the fluid system of the spacecraft. In 
general, the magnitude range for g-jitter accelerations is 
10 -5 to 10 -9 g 0 . 

Figure 6A shows the time-sequence evolution of the 
three-dimensional dynamical behaviour of interface oscil- 
lations driven by the combined gravity gradient and g-jitter 
accelerations and dominated by the gravity gradient force. 
The gravity gradient acceleration given by Equation (14) 
indicates that there are greater positive components of the 
acceleration longitudinal to the direction between the space- 
craft mass centre and the centre of the Earth, and smaller 
negative components of acceleration transverse to this 
direction. As indicated in Equation (15), angle </% varies 
with time. This phenomenon shows that the gravity gradient 
acceleration exerted on the fluid system of the spacecraft 
is equivalent to a combination of the time-dependent twist- 
ing force with turn-around direction and the torsional 
moment acting on the dewar when the spacecraft is orbiting 
around the Earth. Figure 6A shows that the deformation of 
the bubble is created by an asymmetric torsional moment 
with a twisting force. 

The evolution of the capillary-effect-govemed interface 
oscillations at various cross-sections driven by a gravity 
gradient-dominated acceleration is also examined. 
Figure 6B shows the time sequence of the capillary-effect- 
govemed liquid-vapour interface profiles, driven by the 
same gravity gradient-dominated acceleration coupled with 
orbital dynamics as that shown in Figure 6A, in the vertical 
r-z plane at 0 = 0° and 180°. This indicates that the dough- 
nut-shaped bubble (helium vapour) configurations change 
from axial symmetric to asymmetric profiles in both planes 
as i// E ; varies from 0° to 360°. 

The r-z cross-section of the doughnut profiles shows 
oscillations of two bubbles in Figure 6B . Here, the one-up 
one-down and one-down one-up oscillations indicate a very 
important characteristic of the gravity gradient-dominated 
acceleration which produces a combination of a time- 
dependent equivalent torsional moment and a twisting force 
acting on the fluid system of the spacecraft when it is orbit- 
ing around the Earth. 

In addition to the 0.1 Hz low-frequency gravity gradient- 
dominated acceleration, 1 .0 Hz medium-frequency and 
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(B) Interface Oscillations in r-z plane at 0 = 0° and 180° 
(a) t = 0 s (b) t = 534 s (c) l = 884 s 
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Figure 6 Three-dimensional bubble deformations in response to coupling between gravity gradient-dominated accelerations and 
orbital dynamics 


10 Hz high-frequency ^-jitter acceleration-driven liquid- 
vapour interface oscillations have also been investigated. 
Characteristics of these three ranges of ^-jitter frequencies 
were well documented in earlier studies 4 - 5 . Results show 
that lower-frequency g-jitter accelerations contribute more 
in the driving of asymmetric profiles of the interface than 
the higher-frequency g-jitter accelerations. These results are 
in agreement with the earlier studies 45 . 

Figure 7 A shows sloshing due to the time evolution of 
fluid mass centre fluctuations (in the rotational frame) 
driven by gravity gradient-dominated acceleration coupling 
with orbital dynamics. The values of the maximum and 
minimum locations of the fluid mass centre fluctuations are 
Max (* LCf y LC , z LC ) = (0.74, 0.72, 1.35) cm and Min (jr LC , 
Vlo Zlc) = (-0.73, -0.67, -0.01) cm, respectively, while 
the fluctuating values of the maximum fluid mass centre 
disturbances are Max (A x LC , Ay LC , Az LC ) = (1.4, 1.4, 1.3) 
cm. This shows that Ax lx ~ Ay LC > A z LC . Results of fluid 
mass centre fluctuations driven by gravity gradient-domi- 
nated acceleration coupling with orbital dynamics can be 
summarized as follows: (a) fluid mass centre oscillations 
along the jc- and y-directions basically follow the trends of 
the time-dependent evolution of gravity gradient acceler- 
ations along the jc- and y-axes in a tidal mode manner 19 23 
of coupling with orbital dynamics, as shown in Figures 58- 
a and 5B-b\ (b) there are fluid mass centre fluctuations with 
a smaller amplitude along the z-axis in response to the trend 
of the time-dependent evolution of the gravity gradient 
acceleration along the z-axis in a similar manner; (c) coup- 


ling between slosh and orbital dynamics also results in an 
asymmetric tendency of a tidal mode nature 19 23 of the fluid 
mass centre fluctuations while the spacecraft is orbiting 
around the Earth; (d) there is a time lag of 150 s between 
the moment when the gravity gradient-dominated acceler- 
ations act on the spacecraft and the time when fluid mass 
centre fluctuations reach their maximum - in other words, 
the maximum of the fluid mass centre fluctuations was not 
reached instantly and did not change its direction at the 
moment the gravity gradient force changed direction; (e) 
Figure 78 shows the time history of the liquid mass centre 
fluctuations. Large-amplitude fluctuations correspond to the 
directional change of the gravity gradient with respect to 
the azimuth angle variation, while the small ripples 
imposed on the oscillations are contributed by the slosh 
dynamics. 

Spacecraft system mass centre fluctuations (in a non- 
rotational frame) in response to the coupling between slosh 
dynamics driven by a gravity gradient-dominated force and 
orbital dynamics deviating from normal operation have 
been studied, and are illustrated in Figure 7C. The values 
of the maximum and minimum locations of spacecraft mass 
centre fluctuations are Max (X sc , T sc , Z sc ) = (0.35, 0.29, 
0.60) cm and Min (X se , K sc , Z sc ) = (-0.34, -0.32, -0.07) 
cm, respectively, while the fluctuating values of the 
maximum spacecraft mass centre disturbances are Max 
(AX SC , A T sc , AZ SC ) = (0.69, 0.61, 0.67) cm. This shows that 
AX SC ~ A T sc ~ AZ SC in a non-rotational frame. The results 
are very similar to those concluded from the fluid mass 
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Figure 7 Time evolution of fluid and spacecraft mass centre fluctuations driven by gravity gradient-dominated acceleration and 
orbital dynamics. (A) Fluid mass centre in (x, y, z) components (rotational frame), (B) fluid mass centre (absolute value), (C) spacecraft 
mass centre in (X, Y, Z) components (non-rotational frame) and (D) spacecraft mass centre (absolute value) 


centre fluctuations obtained in the previous section, except 
that (a) the former case is in a non-rotational frame and 
the latter case is in rotational frame and (b) the former 
case has smaller fluctuations due to the effect of averaging 
between the dry mass of the spacecraft and the fluid mass 
fluctuations. Figure 7D shows the time-dependent variation 
in the absolute value of the spacecraft system mass centre 
fluctuations. Comparison between Figures 7B and 7D 
shows that the spacecraft mass centre fluctuations are 
mainly dominated by the fluid system mass centre fluctu- 
ations. Failure to consider the effect of slosh dynamics act- 
ing on the spacecraft may lead to the wrong results for the 
case of a gravity gradient-dominated acceleration acting on 
the orbiting spacecraft. 

Coupling between slosh dynamics driven by 
combined gravity gradient and g -jitter 
accelerations dominated by g-jitter forces and 
spacecraft dynamics 

In this case, combined gravity gradient and g-jitter with 
an amplitude of 10“ 6 g 0 for g shown in Equation (17) is 
considered. The gravity gradient acceleration in this case 
is show n in Figure 5B , while the g-jitter acceleration with 
background amplitude of 10 8 g 0 , shown in Figure 5C, is 
replaced by 10“* g n for the present case. Then the combined 
gravity gradient (of the order of 10 7 gj and g-jitter acceler- 
ations acting on the spacecraft correspond to the case of 
orbital accelerations dominated by g-jitter forces acting on 
the fluid system of the spacecraft. As indicated earlier, g- 
jitter accelerations are in the range 10“* to I0 _y g u . 

Figure HA shows the time sequence evolution of the 
three-dimensional dynamic behaviour of interface oscil- 
lations driven by combined gravity gradient and g-jitter 


accelerations dominated by a g-jitter force of 0.1 Hz low- 
jitter frequency. The g-jitter acceleration, shown in Equa- 
tions (16) and (17) and illustrated in Figure 5C, indicates 
that there is a sinusoidal oscillation longitudinal to the 
direction from the spacecraft mass centre to the centre of 
the Earth (parallel to unit vector f c shown in Figure 5 A ). 
As indicated in Equation (15), the azimuth angle t/v varies 
with time. These phenomena show that the g-jitter acceler- 
ation exerted on the fluid system of the spacecraft is equiv- 
alent to the oscillatory forces that push the bubble in the 
combined directions of down-up (see ^-component of g- 
jitter acceleration shown in Figure 5C) and left-right (see 
.v- and y-components of g-jitter accelerations shown in 
Figure 5C as the bubble rotates with respect to the space- 
craft spinning axis. Time variations in the jc- and y-compo- 
nents of the g-jitter accelerations, shown in Figure 5C, 
characterize the left-right oscillations of bubbles shown in 
Figures HA and HB . Figure HB shows the time sequence of 
the liquid- vapour interface profiles in the vertical r-z plane 
at 6 = 0° and 1 80°. Among the various frequencies of the 
g-jitter-dominated accelerations considered, it is shown that 
lower-frequency g-jitter accelerations contribute more to 
the driving of asymmetric profiles of the interface than 
higher frequency g-jitter accelerations 32 3 *. 

Figure 9A shows sloshing due to the time evolution of 
the fluid mass centre fluctuations (in a rotational frame) 
driven by g-jitter-dominated acceleration coupling with 
orbital dynamics. The values of the maximum and mini- 
mum locations of the fluid mass centre fluctuations are Max 
('V, x’> Vlo Zix ) = (1.19, 1.21, 3.53) cm and Min (jc lo y LC , 
Z\ .c) = (-1.20, -1.03, -1.84) cm, respectively, while the 
fluctuating values of the maximum fluid mass centre dis- 
turbances are Max (Ax, r . Ay, c , A z lr ) = (2.38, 2.24, 5.37) 
cm. This shows that Z\ ,c > x LC >Vix ■ Results of fluid mass 
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(B) Interface Oscillations in r-z plane at 0 = 0° and 1 80° 
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Figure 8 Three-dimensional bubble deformations in response to coupling between g-jitter-dominated accelerations and orbital 
dynamics 


centre fluctuations driven by g-jitter-dominated acceleration 
coupling with orbital dynamics can be summarized as fol- 
lows: (a) fluid mass centre oscillations along the x- and in- 
directions basically follow the trends of the time-dependent 
evolution of g-jitter accelerations in the x- and y-axis coup- 
ling with orbital dynamics shown in Figures 5C-a and 5C- 
b\ (b) fluid mass centre fluctuations have greater amplitudes 
along the z-axis in response to the trend of the time-depen- 
dent evolution of a g-jitter acceleration along the z-axis in 
which one half of the orbit period is in the positive and the 
other in the negative direction - this result is exactly 
reflected by the fluid mass centre fluctuations along the z- 
axis shown in Figure 9A with small ripples imposed by 
slosh dynamics; (c) coupling between slosh and orbital 
dynamics also results in the asymmetric nature of the fluid 
mass centre fluctuations while the spacecraft is orbiting 
around the Earth; (d) there is a time lag of 150 s between 
the moment when the g-jitter-dominated accelerations act 
on the spacecraft and the time when maximum fluid mass 
centre fluctuations are reached - in other words, the fluid 
mass centre fluctuations react rather slowly and do not 
instantly change direction when the g-jitter direction 
changes; (e) Figure 9B shows the time-dependent variation 
of the absolute value of fluid mass centre fluctuations 
resulting from coupling between the slosh and orbital 
dynamics. It shows time-dependent fluctuations with large 
amplitudes corresponding to directional changes in g-jitter 
with respect to azimuth angle, while the small ripple 
imposed on the oscillations is contributed by the slosh 
dynamics. 


Spacecraft system mass centre fluctuations (in a non- 
rotational frame) in response to the coupling between slosh 
dynamics driven by a g-jitter-dominated force and orbital 
dynamics deviating from normal operation have been stud- 
ied, and are illustrated in Figure 9C. The values of the 
maximum and minimum locations of the spacecraft mass 
centre fluctuations are Max (X sc , K sc , Z sc ) = (0.29, 0.43, 1 .6) 
cm and Min (X sc , K sc , Z sc ) = (-0.21, -0.21, -0.83) cm, 
respectively, while the fluctuating values of the maximum 
spacecraft mass centre disturbances are Max (AX SC . A T sc , 
AZ SC ) = (0.5, 0.7, 2.5) cm. This shows that AZ SC > A Y sc > 
AX SC in a non-rotational frame. Results of spacecraft mass 
centre fluctuations driven by coupling of a g-jitter-domi- 
nated force and orbital dynamics are very similar to those 
driven by gravity-gradient dominated force and spacecraft 
dynamics. Failure to consider the effect of slosh dynamics 
acting on the spacecraft may lead to the wrong results for 
the case of a g-jitter-dominated acceleration acting on the 
orbiting spacecraft. 

Characteristics of various orbital accelerations con- 
sidered in this study, in addition to the results obtained from 
the previous section for He II slosh dynamics driven by 
various accelerations coupled with orbital dynamics, are 
summarized and illustrated in Table I . 


Conclusions 

The coupling of motions in a spinning spacecraft dewar 
activated by (a) lateral impulse, (b) gravity gradient-domi- 
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Figure 9 Time evolution of fluid and spacecraft mass centre fluctuations driven by ^-jitter-dominated accelerations and orbital 
dynamics. (A) Fluid mass centre in (x, y, z) components {rotational frame), (B) fluid mass centre (absolute value), (C) spacecraft 
mass centre in (X, Y, Z) components (non-rotational frame) and (D) spacecraft mass centre {absolute value) 


nated and (c) ^-jitter-dominated accelerations with space- 
craft orbital dynamics has been investigated. It is shown 
that coupling between He II slosh dynamics and orbital 
dynamics can cause large-amplitude fluctuations to act on 
a spacecraft, resulting in deviations from normal operation. 

For the purpose of this study, the slosh dynamics are 
based on a rotational frame, while the orbital dynamics are 
associated with a non-rotational frame. Computation of He 
II slosh dynamics and the fluctuations of the spacecraft 
mass centre, driven by the three types of force, indicate 
that the major contribution of the orbital dynamics is driven 
by those forces coupling w ith slosh dynamics. Among these 
three environmental and impulsive forces acting on the 
spacecraft, both the gravity-gradient and g-jitter forces are 
in the form of a continuous spectrum while the impulse is 
in the form of a spotty 5-function. All these forces are cap- 
able of exciting slosh dynamics coupling with orbital 
dynamics. In fact, failure to consider the effect of slosh 
dynamics acting on the spacecraft may lead to the wrong 
results. The asymmetric orientation of the helium II liquid 
sloshing in the tank driven by the three forces considered 
may also create control and pointing problems causing 
deviation from normal spacecraft operation 5 i9 . 
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